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ABSTRACT 

We develop a model for the metal enrichment of the intracluster medium (ICM) that 
combines a cosmological non-radiative hydrodynamical TV-Body /SPH simulation of a 
cluster of galaxies, and a semi-analytic model of galaxy formation. The novel feature 
of our hybrid model is that the chemical properties of the diffuse gas in the underlying 
simulation are dynamically and consistently generated from stars in the galaxies. We 
follow the production of several chemical elements, provided by low- and intermediate- 
mass stars, core collapse and type la supernovae. We analyse the spatial distribution 
of metals in the ICM, investigate the way in which the chemical enrichment proceeds, 
and use iron emissivity as a tracer of gas motions. The similar abundance patterns 
developed by O and Fe indicate that both types of SNe pollute the ICM in a simi- 
lar fashion. Their radial abundance profiles are enhanced in the inner 100 h~ x kpc in 
the last Gyr because of the convergence of enriched gas clumps to the cluster centre; 
this increment cannot be explained by the metal ejection of cluster galaxies which is 
quite low at the present epoch. Our results support a scenario in which part of the 
central intracluster gas comes from gas clumps that, in the redshift range of z ~ 0.2 
to ~ 0.5, have been enriched to solar values and are at large distances from the cluster 
centre (from ~ 1 to ~ 6/z -1 Mpc) moving at very high velocities (from ~ 1300 to 
~ 2500 km s -1 ). The turbulent gas motions within the cluster, originated in the inho- 
mogeneous gas infall during the cluster assembly, are manifested in emission- weighted 
velocity maps as gradients that can be as large as ~ 1000 kms -1 over distances of a 
few hundred kpc. Gradients of this magnitude are also seen in velocity distributions 
along sight lines through the cluster centre. Doppler shifting and broadening suffered 
by the Fe K a 6.7 keV emission line along such sightlines could be used to probe these 
gas large-scale motions when they are produced within an area characterised by high 
iron line emissivity. 
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1 INTRODUCTION 

Clusters of galaxies are the largest virialised objects in the 
Universe. Within the hierarchical structure formation sce- 
nario of the cold dark matter (CDM) model, they arise from 
the rare highest peaks of the initial density perturbation 
field. During cluster assembly, the dissipative gas component 
is shock heated by the potential well of the cluster, reaching 
very high temperatures (T ~ 10 7-8 K). This hot and diffuse 
gas radiates energy through thermal bremsstrahlung in the 
X-ray band of the electromagnetic spectrum. The gravita- 
tionally dominating collisionless dark matter component on 
the other hand shows considerable substructure, as revealed 
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by h igh resolution simu lations that follow cluster assembly 
fe.g. lSpringel et aLll200ll and references therein). These sub- 
structures are the surviving remnants of dark matter haloes 
that have fallen in at earlier times, and plausibly mark the 
location of the cluster galaxies. 

X-ray observations have first revealed the presence of 
hot gas trapped in the potential well of galaxy clusters, 
and provide a wealth of information about the dynamical, 
thermal and chemical properties of clusters. A strong fea- 
ture in their spectra is the presence of line emission pro- 
duced by highly ionised iron, mainly, the H- and He-like 
iron lines at 6.9 and 6.7 keV. The content of iron in the in- 
tracluster m edium (ICM) is app roximately one third of the 
solar value (jFukazawa et alJ Il998l : Ettori, Allen & Fabian 
2001), suggesting that some of its hot gas must have orig- 
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inated in the galaxies that reside in the cluster. The pro- 
cesses considered for the supply of this enriche d gas include 
galactic winds driven by supernovae explosi ons (IWhitell99 
Renzh3 Il997h . ram pressure stripping (Mori & Burk 
2000), early enrichmen t by hypernoyae ass ociated with pop- 
ulatio n type HI stars JLo ewenst ei nl 1200 lh. and intracluster 
stars (|Zaritsk Yl Gonzalez fc Zabludofl||2004|) . 

Spatially resolved maps of the temperature and 
metal abundance distributions in the intracluster medium 
have been obtained by X-ray surveys carried out with 
ROSAT, ASCA, BeppoSAX (Finoguenov, David & 
Ponman 2000; Fino g uenoy , Arnaud & David 2001; 
iDe Grandi fc Molendil 1200111 . and more recently with 
XMM - Newton and Chandr a iGastald ello Molendil 
20021: iFinoguenov et ail 120021: ISanders fc Fabianl 120021: 

" 120031 : " 



Matsushr^t^hT^gu^nov BohringerT 2003; ^miuraetliLl 
2004). As a result, abundance profiles of several chemical 
species with different origins have been obtained, namely 
for O and Si generated by core collapse supernovae (SNe 
CC) and Fe, mainly produced by type la supernovae (SNe 
la). These profiles reflect physical mixing processes in the 
baryonic component due to gas cooling, star formation, 
metal production, and energetic and chemical feedback, all 
coupled to the hierarchical build up of the galaxies and the 
clusters. These observations can therefore be very valuable 
for constraining models of galaxy formation and chemical 
enrichment of the ICM, allowing to evaluate the relative 
importance of different types of supernovae for metal 
pollution. The spatial distribution of metals obtained from 
models of the ICM chemical enrichment can be analysed by 
constructing X-ray weighted metal maps which trace the 
interactions between the ICM and cluster galaxies; thus, 
when they are compared with observed metal maps, they 
provide information about the mechanisms involved in the 
enrichment pro cess and may reveal t he dynamical state of a 
galaxy cluster (|Kapferer et al.lEo05l : iDomainko et al.l l2*Q05") . 

X-ray observations also provide important information 
about the dynamics of the intracl uster gas through radial ve - 
locity measurements in clusters (|Dupke Sz Bregmanll200 3) . 
allowing in principle further tests of structure formation 
models. Gaseous bulk flows affect the characteristics of X- 
ray spectra of the ICM by Doppler broadening and shifting 
of emission lines. The prominent Fe K a 6.7 keV iron emis- 
sion line present in X-ray spectra of galaxy clusters could 
be used as a tracer of thes e mo tions. This possibility was 
evaluated bv lSunvaev et all (|2003l) . who showed that future 
X-ray missions like CONSTELLATION-X and XEUS, with 
an energy resolution of a few eV, should be able to detect 
several Doppler-shifted components of this emission line in 
the core of the cluster. 

All this wealth of data calls for a theoretical interpre- 
tation of the ICM chemical enrichment in the framework of 
a detailed cluster formation model. The approaches used so 
far to this end include semi-analytic models of galaxy forma- 
tion coupled to N- Body simulatio ns of galaxy clusters (D e 
Lucia, Kauffmann & White 2004: iNagashima et aflbOOSb . 
hydrodynamical simulations of cluster formation which in- 
clude, self-consistently, star formation, S Ne feedback and 
metal enrichment from SNe CC and la (e.g. IValdarninil 2003 : 
iTornatore et alJl2004Tl . and an intermediate approach that 
combines TV-Body and hydrodynamic simulations together 
with a phenomenological galaxy formation model and a pre- 



scription of the effect under study, such as galacti c winds and 
merger- driven starb ursts (iKaofe rer et a lT 2005) and ram- 
pressure stripping ()Schindlere^nTl2005[ IDomainko et all 
2005). In this work, we present a different intermediate ap- 
proach which combines a cosmological non-radiative hydro- 
dynamical TV-Body/SPH simulation of a cluster of galax- 
ies, with a semi-analytic model of galaxy formation. The 
simplicity of the semi-analytic model has the advantage of 
reaching a larger dynamic range than fully self- consistent 
hydro- simulations, at a far smaller computational cost. In 
particular, it allows us to explore more easily the range of 
parameters that characterise appropriate chemical enrich- 
ment models. 

The long cooling time of the bulk of the gas in ri ch clus- 
ters justifies the assump tion of a non-radiative gas (jEvrardl 
1990: Iprenk et al.lll999f) . However, this time-scale becomes 
smaller than the Hubble time in the core of the cluster, 
where the densities are higher. Here is where the semi- 
analytic model plays an important part, taking into account 
radiative cooling, and models for star formation, chemi- 
cal enrichment and energetic feedback from galaxies. The 
semi-analytic model used here is based o n previous works 
fepringel et al . 11200 ltlPe Lucia et alJl2004 DL04 hereafter), 
but was extended with a new chemical implementation that 
tracks the abundance of different species resulting from mass 
loss through stellar winds of low-, intermediate- and high- 
mass stars, and from the explosions of SNe CC and la su- 
pernovae. The principal new feature of our model is how- 
ever the link between semi-analytic model results and the 
chemical enrichment of gas particles in the underlying N- 
Body/SPH simulation. We pollute gas particles with met- 
als ejected from the galaxies, consistent with the modelling 
of the semi-analytic model. These metals are then carried 
around and mixed by the hydrodynamic processes during 
cluster formation. Thus, we are not limited to an analysis 
of mean chemical properties of the intracluster gas, like in 
the original model of DL04. Instead, the enrichment of gas 
particles allows a study of the spatial distribution of metals 
in the ICM, and to use iron emissivity as a tracer of gas 
motions. 

The aim of our model is, on one hand, to understand 
the way in which the chemical enrichment patterns char- 
acterising the intracluster gas develop, and, on the other, 
to detect the imprints of gas bulk motions in the shape of 
X-ray emission lines, thus providing guidance for the inter- 
pretation of future spectroscopic X-ray data. The first part 
of the present study is therefore devoted to an analysis of 
chemical enrichment of the ICM, focusing on the spatial dis- 
tribution of different chemical elements. In a second part we 
analyse the link between the occurrence of multiple com- 
ponents in the Fe K a 6.7 keV emission line and the range 
of velocities of the gas that generates this emission. This 
is directly relevant to the question whether metal emission- 
line-weighted velocity maps of the cluster could be used to 
infer the spatial coherence of these bulk motions. Gathering 
all the information available in the form of projected maps 
and detailed spectra along lines of sight through the ICM, 
we are able to extend the information in 2-D provided by 
observed projected distributions of gas temperature, abun- 
dances and surface brightness, into a 3-D picture of the ICM 
properties. 

This paper is organised as follows. Section |2] describes 
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our hybrid model used to study the ICM chemical enrich- 
ment, summarises the properties of the hydrodynamical sim- 
ulation used and presents the main features of the semi- 
analytic model. Section contains the details of the chem- 
ical model implemented in the semi-analytic model. It also 
describes the way in which metals released by stellar mass 
losses and supernovae explosions are spread among gas par- 
ticles of the Af-Body/SPH simulation. In Section^] we fix the 
parameters that characterise our version of the semi-analytic 
model and compare its results with local observations. The 
following three sections present the main results of this work. 
Section [5| gives the radial abundance profiles of different 
chemical elements, comparing them with observations. In 
Section |S| we visualise the thermodynamical and chemical 
properties of the ICM by projecting the corresponding mass- 
weighted or emission- weighted related quantities. These two 
techniques (radial profiles and projected maps) are used to 
analyse the history of the ICM chemical enrichment as well 
as the dynamical evolution of the intracluster gas, as dis- 
cussed in Section |7| Section |H| describes the use of Fe K a 
6.7 keV emission line as a probe of the ICM dynamics. Fi- 
nally, in Sectional we give a summary and discussion of our 
results. 

Throughout this paper we express our model chemical 
abundances in terms of the solar photospheric ones recently 
recalibrated by Asplund, Grevesse & Sauval (2005). When 
necessary, we refer the chemical abundances obtained in dif- 
ferent observational works to these standard solar values. 



2 HYBRID MODEL OF THE ICM CHEMICAL 
ENRICHMENT 

Our hybrid model for studying the chemical enrichment 
of the intracluster medium consists of a combination of a 
non-radiative TV-Body/SPH simulation of a galaxy cluster 
and a semi-analytic model of galaxy formation. Dark mat- 
ter haloes and substructures that emerge in the simulation 
are tracked by the semi- an alytic code and used to gener- 
ate the galaxy population fcpringel et al.ll200lTl . The key 
additional component of our model is that we also use the 
diffuse gas component of the hydrodynamical simulation to 
follow the enrichment process. By enriching the gas parti- 
cles locally around galaxies, we can account for the spreading 
and mixing of metals by hydrodynamical processes, thereby 
obtaining a model for the evolution of the spatial distri- 
bution of metals in the ICM. This in turn allows a more 
detailed analysis of the ICM metal enrichment than possi- 
ble with previous models bas ed on semi-analytic techniques 
teauffmann fc Charlotlll 998: DL04). In this section, we pro- 
vide a detailed description of both parts of the modelling. 



2.1 Non-radiative iV-Body/SPH cluster 
simulation 

In the present study, we analyse a non-radiative hydro- 
dynamical simulation of a cluster of gal axies. The cluster 
was s elected from the GIF-ACDM model (|Kauffmann et all 
1999), characterised by cosmological parameters Q o =0.3, 
Q A =0.7 and H G = 100 /ikrns" 1 Mpc -1 , with h = 0.7; it 
has spectral shape F = 0.21, and was cluster- normalised to 



as = 0.9. The second most massive cluster in this parent cos- 
mological simulation (of mass M vir = 8.4xl0 14 /i _1 M©) was 
resimulated with higher resolution in the Lagr angian region 
of th e object and its immediate surroundings (|Springel et alJ 
I2001T) . with a number N hr 1999978 of high resolution 
particles. Within this region, each particle was split into 
dark matter (radm = 1.18 x 10 9 /i _1 M ) and gas (m gas = 
1.82 x 10 8 /i _1 M ) according to Q b = 0.04 and Q D m = 0.26, 
consistent with the baryon density Qbh 2 = 0.02 required by 
Big- Bang nucleosynthesis constraints. However, the identifi- 
cation of dark matter haloes for the semi-analytic model was 
based only on the dark matter particles, with their mass in- 
creased to its original value. The boundary region, where 
mass resolution degrades with increasing distance from 
the cluster, extends to a total diameter of 141.3 h~ x Mpc. 
The simulation was carried out with the paralle l tree N- 
Body/SPH code GADGET (ISoringel et al.1 l200ll : ISpringell 
2005). The starting redshift and the gravitational soften- 
ing length are z s tart = 50 and e = 3.0 kpc, respectively. 
We note that this simulation did not include gas cooling. 
This is a reasonable approximation for cluster simulations 
since the bulk of the halo gas has very long cooling times. 
However, the cooling time becomes shorter than the age of 
the Universe in a central region with radius ~ 100 — 200 kpc, 
from where the bulk of X-ray radiation is emitted (|Sarazinl 




2.2 Semi-analytic model of galaxy formation 

The properties of galaxies in the semi-analytic model are 
determined by the included physical processes. We consider 
the effects of cooling of hot gas due to radiative losses, star 
formation, feedback from supernovae explosions, metal pro- 
duction, and merging of galaxies. Except for the metal pro- 
duction, the parametrization of these physical processes cor - 
responds to the model described by S or in ge l et alJ (|2QQlT) . 
However, we modified the model to track the enrichment 
cycle of metals between the different baryonic components 
of the haloes, i.e., cold gas, hot diffuse gas, and stars, in the 
most consistent way possible within our framework. Our spe- 
cific chemical implementation is similar to the semi-analytic 
model discussed by DL04. In the present work, we refine the 
chemical enrichment prescription of DL04 by tracking the 
mass evolution of different chemical elements as provided 
by three kinds of sources: low- and intermediate-mass stars, 
core collapse supernovae, and type la supernovae. The first 
group of sources yields metals through mass losses and stel- 
lar winds. In the following, we briefly summarise the details 
of the semi- analytic code. In Section |3 we then describe the 
new chemical implementation we introduce in this work. 

Based on the merging trees, the mass of hot gas is calcu- 
lated at the beginning of the evolution between consecutive 
outputs of the simulation. We assume that the hot gas al- 
ways has a distribution that parallels that of the dark matter 
halo, whose virial mass changes from one output to another 
due to the hierarchical growth of structure. Once a fraction 
of hot gas has cooled and star formation and feedback pro- 
cesses are triggered, the mass of hot gas is given by 

M hot = / b M vir - ]T [M« llar + M c ( l], (1) 

i 

where M v i r is the virial mass of the dark matter halo, Mhot is 
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the mass of the hot gas associated with it and M st eiiar and 
M co id are the masses of stars and cold gas of each galaxy 
contained in the halo. The virial mass is given by M v i r = 
100H 2 Rl ir /G. This is the mass enclosed by the virial radius 
i?vir, which is defined as the radius of a sphere with mean 
density 200p C rit centred on the most bound particle of the 
group. Here p C rit is the critical density. The virial velocity is 
given by V v 2 ir = GM vir /R vir . 

The mass of hot gas that cools at each snapshot is given 
by the cooling rate 

dM coo i 2 dr coo i 

-ir =4w -if (2) 

where p g is the density profile of an isothermal sphere that 
has been assumed for the distribution of hot gas within a 
dark matter halo, and r coo i is the cooling radius. The local 
cooling time is defined as the ratio of the specific thermal en- 
ergy content of the gas, and the cooling rate per unit volume 
A(T, Z). The latter depends quite strongly on the metallicity 
Z of the hot gas and the temperature T = 35.9(V v ir/kms _1 ) 
of the ha lo, and is represented by the cooling functions com- 
puted bvlSutherland & Donitd 1199^. 
The star formation rate is given by 

dM* _ aMcoid 
dt ~ - gx ' [6) 



dyn 



where £|* n = 0.1i? v ir/VW is the dynamical time of the 
galaxy, and a is a dimensionless parameter that regulates 
the efficiency of star formation. We adopt the variable star 
formation efficiency introduced by DL04 that depends on the 
properties of the dark halo, a(Vwr) = ao (V v ir/220 kms -1 ) n , 
with ao and n as free parameters. 

Each star formation event generates a stellar mass 
AM, = M*AT/iV s teps, where iVstep = 50 are timesteps of 
equal size used to subdivide the intervals between simulation 
outputs, AT, and to integrate the differential equations that 
describe the changes in mass and metals of each baryonic 
component over these timespans. Each solar mass of stars 
formed leads to a number rice of core collapse supernovae. 
This class of supernovae includes those of type Ib/c and II, 
the former being generated when the progenitor looses its 
hydrogen-rich envelope before the explosion. Ordinary type 
II explosions are the most abundant ones however. The en- 
ergy .Esncc released by each core collapse supernova is as- 
sumed to reheat some of the cold gas of a galaxy, with a 
mass of 

4 ?7cc ^sncc 



AM reh eat = ~€ 



V 2 



■AM,, 



(4) 



where e is a dimensionless parameter that regulates the ef- 
ficiency of the feedback process. 

The fate of the reheated gas is quite uncertain. We 
adopt a model in which the transfer from the cold to the hot 
phase does induce galactic outflows, i.e. the reheated mass is 
kept within its host halo. This model is referred to as reten- 
tion model. The resulting cycle of metal enrichment is con- 
sistent with the way in which metals provided by the semi- 
analytic model are deposited in the gas particles around each 
galaxy in the underlying A^-Body/SPH simulation, described 
in Section ET^l The effect of other feedback schemes (ejection 
and wind models) on the chemical enrichment of the ICM 
is discussed by DL04. They found that these three prescrip- 
tions (including the retention model) predict a similar time 



evolution of the ICM metal pollution, which mainly occurs 
at redshifts larger than 1. As we show in the present study, 
gas dynamical processes driven during the cluster assembly 
occurring at present epochs contribute significantly to de- 
termine the spatial distribution of metals in the ICM, very 
likely erasing any signature of the feedback scheme involved. 

In a hierarchical scenario of structure formation, merg- 
ers of galaxies are a natural consequence of the accretion 
and merger processes of dark matter haloes in which they re- 
side. We directly use the simulation outputs to construct the 
merger histories of dark haloes. To this end we first identify 
dark haloes as virialised particle groups by a friend- of- friend 
(FOF ) algorithm. The SUBFIND algorithm (ISoringel et al.1 
2001) is then applied to these groups in order to find self- 
bound dark matter substructures, and the resulting set of 
gravitationally bound structures is tracked over time to yield 
merging history trees. 

In this subhalo scheme, we distinguish three types of 
galaxies when tracking galaxy formation along the merging 
trees. The largest subhalo in a FOF group hosts the 'central 
galaxy' of the group; its position is given by the most bound 
particle in that subhalo. Central galaxies of other smaller 
subhaloes that are contained in a FOF group are referred to 
as 'halo galaxies'. The (sub) haloes of these galaxies are still 
intact after falling into larger structures. The third group of 
galaxies comprises 'satellite' galaxies. This type of galaxies 
is generated when two subhaloes merge and the galaxy of 
the smaller one becomes a satellite of the remnant subhalo. 
These galaxies are assumed to merge on a dynamical time- 
scale with the halo-gala xy of the new subhalo they reside 
in (jSpringel et al.ll200ll : DL04). In these previous versions 
of the semi-analytic model, satellite galaxy positions were 
given by the most-bound particle identified at the last time 
they were still a halo galaxy. In the present work, we instead 
assume a circular orbit for these satellite galaxies with a ve- 
locity given by the virial velocity of the parent halo and de- 
caying radial distance to the corresponding central galaxy. 
This provides a more robust estimate of the position of satel- 
lite g alaxies, consistent with the dynamical friction formula 
used (Bi nnev &; Tremaind ll987'). 

With respect to the spectro-photometric properties of 
galaxies, we apply the calculations made by DL04, con- 
sidering evolutionary synthesis models that depend on the 
metallicity of the cold ga s from which the stars formed 
(Bru zual &; Charlotl Il993r ) . Our morphological classifica- 
tion of galax ie s is based on the criterion adopted by 
ISpringel et alJ (|200lT) . who used 'shifted' values of the 
Hubble- type T of galaxies to obtain a good morphology den- 
sity relation. Thus, we classify as SO galaxies those with 
< T < 5, and as elliptical and spirals galaxies those with 
lower and higher values, respectively. 



3 CHEMICAL IMPLEMENTATION 

One of the main aims of this work is to reach a better un- 
derstanding of how the chemical enrichment of the ICM 
proceeds by establishing a connection between the spatial 
distribution of several properties of the ICM, such as sur- 
face mass density, mass- weighted temperature distribution, 
X-Ray surface brightness and chemical abundances. 

Recent observations of radial abundance pro- 
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files of different elemen ts (|Fm^£]iei]^\L^_a]J 
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De^ranTn^^tlTl20o5 iTamura e^uTl2004h provide valuable 
constraints on the physical processes involved, especially 
those related to the feedback mechanisms that inject metals 
into the diffuse phase. The available observational data 
calls for a model that can explain the spatial distribution 
of different chemical elements. Thus, it is essential to 
implement a chemical model that takes into account several 
different elements, and discriminates between different 
sources, such as SNe CC and SNe la. 

In our model, information about the spatial distribution 
of the ICM properties can be obtained from the gas par- 
ticles in the hydro-simulation. However, the non-radiative 
simulation used only give the thermodynamic properties of 
the diffuse hot gas. In order to obtain the chemical prop- 
erties of gas particles, we have to establish a link between 
the chemical implementation in the semi-analytic model and 
the metal enrichment of the gas particles. In the following, 
we will first describe in detail how we model the enrichment 
cycle of metals among the different baryonic components in 
the semi-analytic model, and then explain how these quan- 
tities are tied to the gas dynamics in the simulation. 



3.1 Enrichment cycle of chemical elements in the 
semi-analytic model 

In our model, stars can contaminate the cold and hot gas 
because of mass loss during their stellar evolution and metal 
ejection at the end of their lives. The hot gas has primordial 
abundances initially (76 per cent of hydrogen and 24 per 
cent of helium) , but becomes chemically enriched as a result 
of the transfer of contaminated cold gas to the hot phase due 
to reheating by supernovae explosions. This chemical enrich- 
ment has a strong influence on the amount of hot gas that 
can cool, since we are using metal dependent cooling rates. 
This process in turn influences the star formation activity 
which is ultimately responsible for the chemical pollution. 

Our semi-analytic model considers mass losses by stars 
in different mass ranges, from low and intermediate stars to 
quasi-massive and massive stars, taking into account stellar 
lifetimes. Massive stars give raise to SNe CC. Ejecta from 
SNe la are also included. The cold gas component of each 
galaxy becomes gradually more chemically contaminated as 
star formation proceeds and, consequently, new stars that 
are formed are progressively more metal rich. This fact calls 
for the use of metallicity dependent stellar yields. However, 
variations of yields with stellar metallicities are very small, 
so we only consider the mass dependence of stellar yields. 
We adopt yields of different chemical elements estimated for 
stars with solar heavy element abundance. 

Thus, for each stellar mass AM* formed, we determine 
the fraction of mass of stars contained in a given mass range 
by assuming an Initial Mass Function (IMF) <3>(M). The 
mass of chemical element j ejected by stars with masses in 
the range centred in mfe, with lower and upper limits m l k 



and m k , is given by 

AM 3 . = [R k X 3 + Y 3 } AM*, 



(5) 



where X 3 is the solar abundance of element j, and R k and 
Y k 3 are the recycled fraction and the yield of element j for 



the k-th mass range, respectively. The recycled fraction is 
given by 



Rk 



h 

-I 



$(M) r k dM, 



with r k defined as the difference between the mass m k of 
the star at birth and its remnant mass after mass loss due 
to stellar winds and/or supernova explosions. The yield of 
chemical element j newly formed is 



Y 3 
1 k 



$(M) p 3 k dM, 



where p k is the stellar yield of a chemical element j produced 
by a star with mass in the range centred around m k . The 
variable j may refer to H, He or chemical species with atomic 
number larger than 2. 

The mass of species j ejected by each SN la is labelled 
?7ij a . In this case, we are dealing with the total mass ejected 
because this type of supernovae leaves no remnants, thus 



AM 4ia )fc = Vi, k ml AM,, 



(6) 



where r]i ak is the number of SNe la per stellar mass formed 
in the mass interval k, which depend on the model adopted 
for this type of SNe. Note that SNe la do not contribute to 
H or He, hence M 3 . fl , is zero when referring to these two 

' ej(Ia) k to 

chemical elements. 

The transport of different chemical species between hot 
gas, cold gas, and stars, can be expressed as variation of 
the mass of the chemical element j present in each baryonic 
component, namely 



AM c J old 



= -AM coo i A J bot + AM reh eat A J cold , 



+AM cool A{ ot - AM* ^ cold(SF) 
+AMij + AM| (Ia) - AM reh eat At 



+ AM *^co.d(SF) 



RA^AM*, 



(7) 



(8) 



(9) 



where A 3 B — Mg / Mb , is the abundance of different chemical 
elements in each baryonic component, which involves the 
total mass of the baryonic component Mb and the mass of 
species j contained in it, Mg. The suffix B may refer to hot 
gas, cold gas or stars; ^oid(SF) denotes the abundance of the 
element j in the cold gas at the time of birth of AM*. The 
above set of equations takes into account the accumulated 
contribution of different mass ranges that affect the baryonic 
components at a given time, as a result of the combination 
of the star formation rate of each galaxy and the return 
time-scale of the ejecta from all sources considered. 



3.2 Injection of metals in the diffuse gas 

We describe now the procedure used to distribute the chem- 
ical elements generated by the galaxies in the semi-analytic 
model among the gas particles of the AT-Body/SPH simula- 
tion. For each snapshot, we identify the gas particles con- 
tained within spheres of radius 100 h~ x kpc centred on each 
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galaxy. Once the set of iVg as gas particles around a given 
galaxy is found, we apply our chemical enrichment model 
to them, based on the semi-analytic computation of galaxy 
formation. The mass of chemical element j transfered to the 
hot phase is evenly distributed among the set of N g 

as gaS 

particles as 

Am J gas = +-3— AM reheat Ai oM . (10) 

■L vgas 

This prescription arises from equation Q, which involves 
the hot gas component of the semi-analytic model. Note 
that the effect of gas cooling on gas particles is suppressed in 
our model because of the use of a non-radiative simulation. 
This process should however not change the mass fraction 
of different species, since cooling would reduce the mass in 
them proportionally to their abundances. 

The increment of the masses of chemical elements in a 
gas particle in principle modifies its total mass. However, the 
mass of gas particles does not change in the iV-Body/SPH 
simulation. Since we are interested in the resulting chemical 
abundances of gas particles after this enrichment process, 
we redefine the masses of species j contained in gas parti- 
cles such that their sum gives the original mass of the gas 
particle m ga s, as given in the simulation, and the mass frac- 
tion of each chemical element keeps the value reached after 
the contribution given by Equation (|10f) has been added. 

The chemical elements injected into gas particles are 
redistributed among neighbouring particles simulating a 
small-scale mixing process. We apply a smoothing technique 
to the chemical properties of the gas particles within the 
virial radius of the cluster, considering 32 neighbours around 
each of them. The radius of the sphere containing these 
neighbouring particles is less than ^ 100/i _1 kpc, therefore 
the large-scale distribution of metals in the ICM is not af- 
fected by this procedure. After smoothing, the mass of chem- 
ical species j of particle i is given by 

32 

m gas,i = ^2 m i a s,i/p^s,iW(rii,hi), (11) 
1=1 

where m J gas l and p gas ,i are the mass of species j of the l- 
th neighbour and its density, respectively. W(ru,hi) is the 
kernel used, ru is the relative distance between particles, 
and hi is the smoothing length of particle i for which the 
smoothed quantity is being estimated. The smoothing kernel 
is the one used for SPH calculations in the AT-Body/SPH 
simulation ([Springe!. Yoshida Sz Whitel200lh , that has been 
obtained from a 5-spline. This pocedure is applied at every 
output of the simulation. 



4 CHARACTERISTICS OF THE 
SEMI-ANALYTIC MODEL 

Our semi-analytic model is characterised by several param- 
eters that regulate the way in which gas cooling, star forma- 
tion, supernovae feedback and galaxy mergers proceed. A 
small number of them are considered free parameters and 
their appropriate values are determined by observational 
constraints. These are the star formation efficiency a and 
the feedback efficiency e, which regulate the star formation 
rate and the core collapse supernovae rate, respectively, and 



a parameter characterising the contribution of supernovae 
type la. This section describes the fixed parameters adopted 
and the conditions we impose to determine the free ones. 
We then compare some properties of our model with obser- 
vational results. 



4.1 Fixed parameters of the model 

Our chemical implementation involves an IMF and stellar 
yields of different species. We do not analyse the impact 
of changing these parameters on the chemical enrichment 
of the ICM; they are considered fixed as well as the baryon 
fraction For the later, we adopt a value of 0.13, which was 
chosen to match that used in the Af-Body/SPH simulation. 

We consider a Salpeter IMF, normalised by fixing the 
fraction £ of stars with masses larger than 1M©, which are 
the major contributors to the chemical enrichment (Porti- 
nari, Chiosi &; Bressan 1998). Thus, we obtain the condition 

M u 

J $(M)MdM = £, (12) 

1M 

with the upper limit M u = 1OOM . The lower limit M\ is de- 
termined by the value of £ adopted, such that the integration 
between M\ and M u equals unity. We adopt £ = 0.5, a high 
enough value to reach the right metallicity in all baryonic 
components. 

For low- and intermediate-mass stars (O.8M < M < 
5 - 8M ), we adopt the stellar yields of IMarigol ll2QQlK . 
and for quasi- massive (5M < M < 8M ) and massive 
st ars (8M < M < 120 M ), we use yields from models 
of IPortinari et alJ (|1998fl . The first two ranges of masses 
eject mainly 4 He, 12 C, 14 N, and possibly 16 O, through stel- 
lar winds. The last two stellar mass ranges contribute to the 
chemical enrichment of the interstellar medium by stellar 
winds, and supernovae explosions triggered by core collapse. 
SNe CC are the main contributors of a— elements. In our cal- 
culations, we follow the production of 8 chemical elements 
(H, 4 He, 12 C, 14 N, 16 0, 24 Mg, 28 Si, 56 Fe) generated by stars 
with masses distributed in 27 mass ranges, from 0.8 M to 
100 M . The mass of the rest of the elements produced is 
stored in a separate variable. 

Since we consider yields from stars with solar metal- 
licity, and taking into account that the solar heavy ele- 
ment abundance has been recalibrated to Z = 0.0122 
(|Asplund et alJ 120051) . we estimate average yiel ds from 
those given for met allicitie s 0.02 and 0.008 by IMarigol 
(|200lT) and IPortinari et alJ (|l998F ). Taking into account 
that the combination of standard IMF and yields do not 
satisfy observational constraints, such as abundances in 
solar neighbourhood stars, metallicity-mass relation and 
ICM metallicity (Moretti, Portinari & Chiosi 2003), and 
considering the uncertainties that a ffect the stellar yields 
(|Gibson. Loewenstein &; Mushotzkvlll997f) . Mg and O are 
increased by a factor of 4 and 1.5, respectively. When all 
mass ranges are considered, we get a net yield of 0.043 
(~ 3.5 Zq), and a corresponding recycled fraction of 0.39. 

Integrating the normalised IMF between 8M and 
100 Mq, which is the mass range of stars leading to SNe CC, 
gives r]cc = 0.009. The energy ejected by each SNe CC is 
set to 1.2 x 10 51 , following the results bv lWooslev &; Weaver! 
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(1995), on which the nucleosynthesis calculations of total 
ejecta of massive stars made bv IPortinari et ah! |l998) are 
based. 

We include the contribution of SNe la which are char- 
acterised by high iron production (~ 0.6 M ) and long delay 
time from the formation of the progenitor to the supernova 
explosion. The currently accepted scenarios for the occur- 
rence of SNe la are, on one hand, carbon deflagration in C-0 
white dwarfs in binary systems (single degenerate systems), 
and, on the other, mergers of two white dwarfs (double de- 
generate systems). The former would account for ~ 80 per 
cent of the events, with mean time delays of 0.5—3 Gyr, while 



(ii) e -- 



= 0.1, 
0.2. 



lower time delays of ~ 0.3 G 


vr (iTutukov &; Yungelson 




Yoshii. Tsuiimoto & Nomot< 




Ruiz-Laouente & Canal 


1998: Hachisu. Kato & Nomotnl 1199* iGrecnrio 


2005). We 



adopt the single degenerate m odel to estimate the SN e 
la rate, following the scheme of iGreggio &; Renzinil Jl983), 
where type la supernovae originate in binary systems 
whose components have masses between 0.8 and 8M . 
Calculations are based on t he formalism described in 
iLia. Portinari Carrarcl h()02). but assuming that SNe CC 
originate from single stars with masses larger than 8M . 
For the chemical contribution of SNe la, we adopt the nu- 
cleosynthesis p r escrip tions from the updated model W7 by 

llwamoto et aljflaai)- 

We use the st ellar lifetime given by 
IPadovani Sz Matteuccil (|l993h to model the return time- 
scale of the ejecta from all sources considered, being 
specially relevant for the single stars in the low- and 
intermediate-mass range and for the progenitors of SNe 
la. In the last case, the mass range of secondary stars in 
binary systems gives explosion times for SNe la comprised 
between - 2.9 x 10 7 and - 1.4 x 10 10 yrs, with the SNe 
la rate reaching a maximun within ~ 0.1 — 0.7 Gyr for a 
single stellar population. 

This chemical model shares similar characteristics with 
that implemented in th e semi-analytic code developed by 
iNagashima et ahl (|2005l) , being an improvement with respect 
to other semi- analytic models j Kauffman n &; Char lotl l 1998 ; 
ICole et aLlboont ISomerville et al.ll200ll . DL04) which typi- 
cally consider both the fraction of metals ejected by the stars 
formed and the recycled mass that has been locked in stars 
at their birth and re- ejected as free parameters. 



4.2 Parameters governing the SF and SNe CC 

The parameters involved in the star formation and feedback 
processes have a strong influence on the star formation rate. 
The SNe CC rate closely resembles the star formation rate 
density (SFR) because it is dominated by high mass stars 
with short lifetimes (< 3 x 10 7 Gyrs). The value of the index 
n involved in the star formation law adopted (Equation 
is set equal to the value suggested by DL04, n — 2.2. How- 
ever, the values of ao and e used by DL04 for the retention 
model do not constitute the best option for our semi-analytic 
model. The different chemical model implemented in our 
study demands a new determination of the appropriate val- 
ues. We use as observational constraints the Milky Way 
properties, the luminosity function, the T ully- Fisher, color- 
magnitude and mass-meta llicity relations ( Kauff mann et all 
J l200li DL04), obtaining the best fit for 



The agreement of our model results with the observational 
data used to constrain these free parameters is quite similar 
to that showed by the models of DL04, therefore we refer the 
reader to their paper for further discussion on this topic. 

The evolution of SFR in our model is broadly con- 
sistent with the observ ational results recently compiled by 
ISomerville et all (|200lh . with the peak of the SFR however 
shifted to higher redshifts. This behaviour is similar to re- 
sults for the c osmic star formation history in hydrodynamic 
simulations bv lSoringel Hernauistl ([2003F ) (see their figure 
12), although the normalisation is less certain in our model 
here since our estimate of the comoving volume is made un- 
certain by the irregular shape of our high-resolution zone 
within the simulated volume, which is used to estimate SF 
and SNe density rates. 

We note that our simulatio n method is really q u ite dif - 
ferent from that applied by Soringel & Her nauistl (2003). 
The present work is based on a non-radiative sim- 
ulation combined with a semi-analytic model, while 
Sor ingel Hernauistl (|2003f) use a fully dissipative simula- 
tion with a subresolution model for the ISM and feedback ef- 
fects. However, despite the huge differences in methodology, 
the principal qualitative results are pretty similar, which is 
an encouraging consistency. Even th ough our simulation is 
simila r in volume to simulation 'G5' of lSoringel Hernauistl 
(2003), the SFR density obtained from the semi- analytic 
model is higher than in 'G5', specially at higher redshifts, 
being in better agreement with their simulations with higher 
mass resolution (see their figure 10). One possibility to ex- 
plain this is that the semi-analytic model is less affected by 
resolution effects all the way down to the mass resolution 
limit of the haloes, while this is not the case for the full 
hydro- simulations. While a halo with ~ 30 particles will be 
able to yield a good representation of the SFR when the 
semi- analytic model is applied, one has to go substantially 
above the halo detection threshold in the full hydrodynam- 
ical simulations in order to recover a numerically converged 
result for the hydrodynamics. Effectively, this means that 
the dynamic range of the semi-analytic method is higher. 

A convergence test of cluster simulations similar to the 
one used in the present study has been done by Ciardi, 
Stoehr & White (2003). They analysed the redshift evo- 
lution of the total star formation rate, normalised to the 
total baryonic mass, obtained by applying a semi-analytic 
mode l to simulations 'S1'-'S4' performed bv ISoringel et ahl 
(|200 J) . They show that simulation 'S2', with the same 
mass resolution than our hydrodynamical simulation, con- 
verged at z ^ 9 to the higher resolution simulation 'S3' 
(mam = 2.4 x 10 8 /i _1 M©) which already accounts for all 
significant star formation. Therefore, the simulation used in 
our hybrid model is suitable for studying the chemical en- 
richment of the ICM. 



4.3 Parameters governing the SNe la rate 

The contribution of SNe la to the metal production of the 
model is characterised by the fraction A of binaries origi- 
nating this type of supernovae. Apart from the iron content, 
this parameter do not significantly affect the properties of 
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the baryonic components in the model, hence the constraints 
imposed for determining the star formation and feedback ef- 
ficiencies are almost insensitive to the value of A. This value 
has to satisfy other kind of requirements. 

During the last five years, great efforts have been un- 
dertaken to detect SNe over cosmological distances, with the 
aim of studyi ng the evolution of SNe CC and SNe la rates 
with redshift (I Pain et alll99d Cappellaro. Evans fc Turatto 
1999, [Hardin et al J 1200ft IPain et al.1 12002: iMadgwick et~ 
200i iTonrv et alJbOOA iDahlen et al J 12004 IStroker e~ 
2004). We adopt A = 0.1 in order to reproduce this set of 
observations and to achieve the observed iron distribution 
in the ICM, that will be discussed in detail in Section |5] 

The combination of the SF history of the model and the 
scenario chosen for the origin of SNe la naturally leads to 
SNe CC and SNe la rates whose ratio evolves with redshift. 
The ratio between the number of SNe la and SNe CC has 
a fairly constant value of ~ 0.4 to z y Q.7, similar to th e 
behaviour suggested by observations of lDahlen et alJ (|2p04). 
The time delay in SNe la explosions makes this ratio smaller 
at higher redshifts. We find that it decreases to ~ 0.25 at z ~ 
2.5, where SNe la rate reaches its maximum, and continues 
decreasing afterwards, manifesting the minor contribution 
of SNe la at large lookback times. 

4.4 Comparison of model results with 
observations 

Having specified the parameters that define our model, we 
make further comparisons with observations in order to 
check that it produces a proper circulation of metals among 
the different baryonic components. 

An important constraint on cluster chemical models is 
the fraction of heavy elements ejected by both types of SNe 
which end up in stars and in the intergalactic medium. A re- 
cent study of metal enrichment of early type galaxies based 
on the X-ray properties obtained w ith Chandra has been 
done by |Humphrey Buote] $MH&- The sample of galax- 
ies span ro 3 orders of magnitude in X-ray luminosity (lot), 
from group- dominant to low-Lx galaxies. We apply this se- 
lection criterium to the central galaxies in our simulation 
and estimate the iron abundance of the hot gas con- 
tained in the dark haloes in which they reside. The three 
panels in Figure Q] present the dependence of Z?e with lumi- 
nosities L x (which calcula tion is based on the formula pre- 
sented bv IValdarninill20o3) and Lb, and the iron metallicity 
[Fe/H] of their stars; they are compa red with the values re- 
ported bv lHumphrev &: Buotd f|2QQ5h . The general trend of 
the properties of our model galaxies is consistent with that 
showed by their observational results. On one hand, there 
is no evidence of a correlation between Zpe and Lx or Lb- 
On the other, the iron abundance for the gas and stellar 
components of each of our model galaxies closely follow the 
ZFe(gas) = [Fe/H] (stars) line, around which many of the 
observed galaxies are clustered. This fact strongly supports 
the circulation of metals among the different baryonic com- 
ponents achieved by our model. 

The fraction of iron in hot gas and stars that is syn- 
thesized by different type of SNe is indicative of the dif- 
ferent time-scales of star formation and of pollution of the 
hot gas with SNe la products. Adopting standard SNe la 
and CC stellar yields and performing a fitting analogous to 



that of iGastaldello Molendil (l2002h. iHumnhrev Buotd 
(2005[) infer an SNIa iron enrichment fraction of ~ 0.66 in 
the hot gas of their observed galaxies, and ~ 0.35 in their 
stars. The model galaxies considered in Figure Q] are char- 
acterised by similar fractions; from the total amount of iron 
originated in SNe la, ~ 65 per cent is contained in the hot 
gas of these galaxies while the remaining ~ 35 per cent is 
locked in the stars. The corresponding fractions for iron pro- 
duced by SNe CC are ~ 62 and ~ 38 per cent for the hot 
gas and stars, respectively. We see that the iron mass frac- 
tion in stars is a bit higher if we consider iron from SNe 
CC instead of that produced by SNe la. This fact reflects 
the increment of the ratio between the number of SNe la 
and SNe CC to lower redshifts, thus the hot gas continues 
being polluted by SNe la ejecta after most of the stars have 
been formed. However, the bulk of the star formation have 
occurred when the hot gas have already been considerably 
polluted by SNe la products, as it is evident from the small 
difference between these fractions. 

A similar conclusion can be obtained from the fraction 
of iron originated in different sources that end up in the 
baryonic components of our simulated cluster. We have that 
~ 78 per cent of the iron produced by SNe la is contained 
in the hot gas and ~ 21 per cent is locked in the stars, while 
for SNe CC we get ~ 72 and ~ 28 per cent for hot gas and 
stars, respectively. The interpretation about the influence 
of the evolution of SNe rates on these percentages becomes 
clearer from the inspection of Figure and the associated 
discussion, which focalise on the ejected iron mass rate for 
both types of supernovae produced by the cluster galaxies. 

Our simulated cluster has cold baryon fraction / c , c is — 
M s teiiar/(M st eiiar + Micm) ~ 0.081, iron yield y c is,Fe = 
M Fe /M ste ii ar = (M[Sm + Cellar) /Stellar ~ 3.2 times solar 
and iron mass fraction fp e = MF e /M v [ r ~ 5.3 x 10 -5 , where 
M s teiiar is the mass of stars of the cluster galaxies, Micm is 
the mass of the intracluster hot gas, and M^| llar and M^ M 
are the masses of iron contained in these two baryonic com- 
ponents. These values are compared with estimates from 
a combination of near-infrared properties of galaxies within 
clusters for which X-ray imaging data is available (Lin, Mohr 
& Standford 2003). The properties of galaxies ar e obtained 
from the Two Micron All Sky Survey f2MASS). ILin et al.1 
(2003) present several relations for low redshift clusters 
(0.016 < z < 0.09), spanning a range of X-ray emission- 
weighted temperatures of 2.1 keV < /cbLx < 9.1 keV (with 
&b the Boltzmann constant) that corresponds to about an 
order of magnitude in cluster virial mass (O.8-9xlO 14 M ). 
The simulated cluster has a virial temperature of Ajb T ~ 
7.2 keV. The approxima te values o f the qu antities mentioned 
previously, obtained bv iLin et alJ (J2003) for a cluster with 
this temperature are / c , c is ~ 0.09, Yds,Fe ~ 3.5 times solar 
and /Fe ~9x 10 -5 . Model values are lower than these ones 
but are within the e rror bars of the observations (see figure 
9 of lLin et aliboofl . 

The number abundance of iron relative to hydrogen for 
the ICM respect to the solar value is [Fe/H]icM ~ 0.28, close 
to the lower limit of obse rved ICM abunda nces which range 
from ~ 0.3 to ~ 0.5 fe.g.. lEttori et al .1120011) . The iron mass- 
to- liglr^atiojs^ 0.014 T©, well within the range obtained 
bv lRenzinil (Il993h for h = 0.7, that is, 0.0085 - 0.017 T©. 

The total 5-band luminosity of the cluster is Lb = 
3.5 x 10 12 L©. The 5-band and ]/-band mass-to- light ratios 
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Figure 1. Hot gas iron abundance, Zp e , as a function of the X-ray luminosity, Lx (left panel), £>-band luminosity of the galaxy, Lb 
(middle panel), and iron abundance of their stars [Fe/H] (right panel). Filled circles represent central galaxies in the simulation selected 
according to the X-ray luminosit i es of their associated hot gas (39.8 < log(Lx) < 43.5), as the range span by the sample of galaxies 
analysed bv iHumohrev &; Buotd J2005). Their observational results are represented by empty squares and triangles, where the later 
indicate lower limit determinations of iron abundances. The orbservational points are consistent with the Zp e (gas) = [Fe/H] (stars) line, 
plotted i n the ri ght panel, w ithin t heir error bars, which are not included in these plots for clarity (we refer the reader to the figures 5 
and 6 of Humphrey BuotdEoOsh. 



are Tb ~ 437 T© and Ty ~ 345 T©, respectively. Observa- 
tional determinations of these quanti ties give 7v ~ 175 — 252 
(ICarlberg et"allll996l: iGirardil boOCb . and j B ~ 200 - 400 
(Kent & Gunn 1983; Girardi 2000). The intra-cluster mass 
to light ratio of our cluster is Micm/^b ~ 4Or , quite sim- 
ilar to the typical o bserved value of no 36F0 indicated by 
iMoretti et aiTfcOOA . 

The analysis of the physical properties of gas particles 
included in the hybrid model provides more detailed features 
of the ICM, regarding the spatial distribution of different 
chemical elements, and the connection between the X-ray 
spectra generated by the diffuse gas with its large-scale mo- 
tions, as we will show in the next sections. 



5 METAL DISTRIBUTION IN THE ICM 

The increase of spectral and spatial resolution of the 
last generation of X-ray satellites, like ASCA and Bep- 
poSAX, and more recently, Chandra and XMM-Newton, has 
provided a great deal of information about cluster cen- 
tres. In particular, it has become possible to obtain de- 
tailed abundance measurements of several chemical elements 
within radial bins, allowing a determination of their spatial 
distribution (IFino guenov et alJl20 0(t iDe G randi & Molendi 
| 200ll: iGastaldello fc Molendi l2002t iDe Grandi et al.1 120031: 
iTamura et alJ I2004T) . Radial distributions of these abun- 
dances support the presence of metallicity gradients in the 
ICM, a feature that has been shown to be common in groups 
and clusters of galaxies. These radial profiles are of great im- 
portance to infer the origin of the different species found in 
the cluster. 

In our approach, the spatial distribution of chemical 
and thermodynamical properties of the ICM is provided by 
the gas particles in the Af-Body/SPH simulation. The hy- 
brid model combines the dynamics of gas particles with the 
chemical information supplied by the semi-analytic model. 
Thus, at each snapshot, gas particles carry the information 



about masses of H, He and heavy elements, with the possi- 
bility of isolating the contribution due to SNe CC and SNe 
la. These data are analysed by estimating radial profiles of 
the chemical elements abundances. 

Radial profiles are constructed by taking into account 
the gas particles lying inside the virial radius of the clus- 
ter. We then divide the volume limited by this radius into 
concentric spherical shells centred on the dominant cluster 
galaxy, each of them containing the same number N^f of 
particles. For each shell, we estimate the mass-weighted and 
X-ray emission- weighted mean abundance relative to hydro- 
gen for several species. The later is estimated using the Fe 
Kq, 6.7 keV emission line, which calculation for each gas par- 
ticle is described in section |H| These quantities are then ex- 
pressed in terms of the solar value. These mean abundances 
are assigned to the mean radius of each shell. 

Figure shows the radial distributions of Fe K a 6.7 
keV emission- line- weighted iron, oxygen and silicon abun- 
dances by number for our model, which are preferred to 
mass- weighted ones in order to compare with observations. 
The radius is normalised with LsW- These results are char- 
acterised by the presence of negative gradients in the ra- 
dial distribution of all the chemical elements considered. 
The emission- weighted values are no 0.1 — 0.15 dex higher 
than mass-weighted abundances, specially in the inner re- 
gion (< O.IR/Rvir). 

There are ambiguous results regarding the de- 
pendence of iron abun dance on cluster temperature. 
Bau mgartner et "all (120051) find decreasing abundances from 
~ 0.7 to no 0.3 (Fe/H)© in the temperature range 
JzbT no 3 — 10 keV, while a roughly con stant value of 
no 0.3 (Fe/H)0 is claimed bv iRenzinil (12003). An interme- 
diate beha viour was detected previo usly bv lFukazawa et all 
f(l998f) and IFinoguenov et alJ (|200lf) . with iron abundances 
of - 0.3 (Fe/H)© for cool clusters and - 0.2 (Fe/H)© for 
hotter ones. Thus, it is better to compare observed and sim- 
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Figure 2. Fe K a 6.7 keV emission-line-weighted abundance pro- 
files of Fe, O and S i relative to H, referred to the solar value 



( Asolund et al. 2005). Total abundances, with the contribution 
of both types of SNe, are given by solid lines; root mean square 
standard deviations are indicated by dotted lines. Separate con- 
tributions of SNe CC and SNe la are shown by dashed and 
dashed-dotted lines, resp ectively. Symbols wi th error bars are 
observational data from iTamura et alJ fl2004T) for hot clusters 
(k^T > 6keV, filled circles) and medium-temperature clusters 
(3 < k^T < 6keV, open diamonds). 



ulated results on cluster abundances within the same tem- 
per ature_£an£e : _^ 

ITamura et all (|200A analysed XMM-Newton observa- 
tions of 19 galaxy clusters. They obtained accurate abun- 
dance distributions from deprojected spectra for clusters in 
different temperature ranges, and for three spherical regions 
with radii between 0-50, 50-200 and 200-500 /i _1 kpc. Our 
simulated cluster has temperatures of the order of k^T ~ 7.5 
keV (see mass-weighted and emission-weighted maps in Fig- 
ures El and 3| . Then, results from our model are compared 
with data for hot clusters (/cbt > 6 keV) and intermediate- 
temperature clusters (3 < k&T < 6keV), represented in 
Figure |21 by filled circles and open diamonds, respectively. 
The later set of clusters is used since they are better deter- 
mined, providing better constraints to t he be haviour of O 
and Si. Following I5e Grandi fc Molendil ( 2001), we estimate 
i?vir for observed data from an ICM temperature-dependent 
formula. 

From Figure |21 we can see that the contribution of SNe 
la (dashed-dotted lines) is necessary to recover the observed 
Fe abundance distribution. The contribution of Si from SNe 
la is not as large as for the Fe production, but helps to reach 
the mean observed values. However, the error bars in the Si 
determinations are much larger than those for Fe and do not 
establish a reliable constraint for the model parameters. This 
is also the case for O; the large observational uncertainties 
in the determination of these chemical abundances for hot 
clusters allow any type of behaviour of the radial profiles, 



while data points for medium temperature clusters are con- 
sistent with a flat profile with a slight increasing trend with 
radius. Note that the model abundances show radial pro- 
files with negative gradients for all the elements considered. 
However, the uncertainties both in observed and simulated 
abundances do not allow us to infer strong constraints from 
the comparison at this point. Only the slopes of iron abun- 
dance profiles are better established by observational data. 
The good agreement between the iron abundance distribu- 
tion in our model and that observed in hot clusters makes 
this model suitable for further analysis. We now carry out 
a deeper inspection of our results in order to gain insight in 
the enrichment history and dynamical evolution of the ICM. 



6 MAPS OF ICM PROPERTIES 

Our hybrid model for investigating the ICM metal enrich- 
ment provides information about the chemical and ther- 
modynamical properties of the diffuse intracluster gas in 
terms of the gas particle properties (positions, velocities, 
density p, and internal energy per unit mass u) of the un- 
derlying hydrodynamic simulation. From these quantities we 
can obtain their temperature as T — /jm p /kB (7 — 1)^, 
where \i — 0.6 is the mean molecular weight, m p the 
proton mass and 7 = 5/3 the adiabatic exponent. We 
can learn about the density, temperature and abundance 
distributions, and velocity structure of the ICM by con- 
structing maps of these properties using planar projections 
of the gas particles within the cluster virial radius. This 
kind of analysis is matched to the recent progress in X- 
ray observations with Chandra and XMM-Newton, which 
provides temperature and abundance maps of cluster cores 
(ISchmidt et al.ll2002l: Zanders Fahianl2002l [Sanders et al.l 
2004; Fukazawa, Kawa no Kawashimal l2004r) . In order to 
compare with observations we will analyse not only mass- 
averaged properties, but also X-ray emission- weighted ones. 
As a result of the high temperature and low density of the 
ICM, the most prominent feature in its X-ray spectrum is 
a blend of emission lines from iron (mainly FeXXV and 
FeXXVI) with photon energies between 6.7 and 7 keV. 

The first step in computing line-emission weighted prop- 
erties of the ICM from our simulations is to estimate the 
intensity of the Fe K a 6.7 keV emission line for each gas 
particle. The atomic processes involved in the generation 
of the emission are quite sensitive to the hydrogen volume 
density n#, temperature T, and chemical a bundances of the 
gas. F ollowing the procedure performed by Furlanetto et 
( 2004), we build a grid where these quantities are varied, and 
we compute the intensity of the emission line for each point 
in the grid. For this purpose we use the radi a tive-collisional 
equil ibrium code CLOUDY (|Ferlandl 1200(1 : iFerland et ail 
2001). Given the thermodynamical conditions of the ICM, 
we can well approximate it as a low density hot plasma in 
collisional equilibrium, without being affected by an ioniz- 
ing backgroud. Therefore, photoionization is discarded in 
the models we generate with CLOUDY. 

The iron line emissivity of the gas particles that make 
up the ICM is obtained by interpolating the emissivities 
assigned to our grid of gas states, where the parameters 
denning the three-dimensional grid in density, tempera- 
ture and iron abundances cover the range of values exhib- 
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Figure 3. Projection of mass-weighted (MW) properties of gas particles contained within a sphere of radius l/i -1 Mpc centred on the 
dominant cluster galaxy: hydrogen number density (top-left panel), temperature (top-right panel), iron abundance by number Fe/H 
referred to the solar value (bottom-left panel) and emissivity of the Fe K a 6.7 keV emission line normalised with the maximum mass- 
weighted emissivity I Q (bottom-right panel). 



ited by the gas inside the virial radius of the cluster at 
z = 0. An important detail of our chemical implementa- 
tion is that we follow the evolution of the abundances of 
several species. Hence, instead of assuming solar metallicity 
for gas particles, or varying this quantity as a global param- 
eter like in other wor ks (e.g., ISutherland &; Dooital ll993: 
iFurlanetto et alJ l2004). our grid of models depends on sets 
of abundances of the chemical elements considered in the 
semi-analytic model. As gas particles become more metal 
rich, the abundances of all elements increase roughly propor- 
tionally to each other. We take advantage of this feature to 
construct the sets of abundances characterised by the mean 
value of iron abundance within a certain range of values. 
Since we only consider the most abundant elements in our 
semi-analytic model, the code CLOUDY automatically gives 
solar v alues (meteroritic abundances of [ G mvesse ^^ Anders] 
(I1989F) with extensions bv lGrevesse fc Noels! (11993^ for the 
rest of the elements. 

The range of values covered by the grid parameters 
are 1 x 10 6 < T < 1 x 10 9 K for gas temperature, and 



3.16 x 10" 6 < riH < 0.1 cm for the hydrogen volume den- 
sity. This last range of values is estimated from the density 
p of gas particles in the A/-Body/SPH simulation and from 
their hydrogen abundance A^ as in our model. In general, 
this quantity ranges from the initial value of 0.76 to ~ 0.70, 
which is reached by the most metal rich particles. Iron abun- 
dances by number range from Fe/H~ 5 x 10 -8 to ~ 5 x 10 -5 , 
which are ~ 1.8 x 10 -3 and ~ 1.8 times the solar value, re- 
spectively. 

We choose density spacings A log(riH/cm _3 ) = 0.25. 
Since emissivities are quite sensitive to temperature varia- 
tions, we adopt temperature spacings of Alog(T/K) = 0.1, 
for 10 7 K < T < 10 8 5 K, and Alog(T/K) = 0.25 otherwise. 
We consider seven bins for the iron abundances with spac- 
ings A log (Fe/H) = 0.5 from log(Fe/H) = -7.75 to -4.25. 
In order to produce the whole grid of models in a systematic 
way, we use a user-friendly package of IDL routines that fa- 
cilitates the use of the code CLOUDY (in version C96B4), 
called MICE (MPE IDL Cloudy Environment), developed by 
Henrik Spoon. 
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Figure 4. Projection of Fe K a 6.7 keV emission line- weighted (EW) properties of gas particles contained within a sphere of radius 
1 Mpc centred on the dominant cluster galaxy: hydrogen number density (left panel), temperature (middle panel) and iron abundance 
by number Fe/H relative to the solar value (right panel). 



Figure |3] shows the projected mass- weighted maps of 
hydrogen density (top- left panel), temperature (top-right 
panel), iron abundance (bottom- left panel) and emissivity of 
the Fe K a 6.7 keV line normalised with the maximum mass- 
weighted emissivity I Q (bottom-right panel) for gas parti- 
cles within a radius of 1 h~ x Mpc from the cluster centre at 
z = 0. From the first three maps we can connect the iron 
abundances with the distributions of gas density and tem- 
perature in the ICM. We can see that those regions with 
high levels of chemical enrichment are associated with cool 
and high density areas; this fact is particularly evident in the 
small substructures in the outskirts of the cluster that have 
been recently accreted. These substructures do not produce 
strong X-ray emission due to their low temperature (fourth 
map), and the line emission map results in rather smooth 
maps, following mainly the coarse features of the density 
distribution. 

Figure0]is equivalent to Figure^ but the Fe K a 6.7 keV 
emission- line- weighted quantities are shown instead. We can 
see that the denser, cooler and more metal rich substructures 
found in the outskirts of the mass- weighted maps of Figure |3| 
do not appear in the corresponding plots of Figure 0] since 
the intensity of the lines emitted from these regions is very 
low (see bottom-right panel of Figure |3| ■ 



7 ENRICHMENT HISTORY AND 

DYNAMICAL EVOLUTION OF THE ICM 

In order to understand the way in which the chemical enrich- 
ment patterns characterising the intracluster gas develop, 
we analyse the temporal evolution of maps of the projected 
chemical properties of the gas and the corresponding radial 
abundance profiles. This allows us to evaluate the relative 
contribution of SNe CC and la to the metal enrichment, as 
well as to investigate the dynamical evolution of the intr- 
acluster gas which plays a fundamental role in the spatial 
distribution of metals. 

For this purpose, we first consider gas particles lying 
within 500 h~ x kpc from the cluster centre at z — and trace 



their properties back in time. Thus, we follow the chemical 
enrichment of small substructures forming at high redshift 
and study how they progressively become part of the cluster 
and contribute to establishing its final distribution of met- 
als. Figures |5] and El show mass- weighted iron abundances by 
number relative to the sun, Fe/H (left column), the corre- 
sponding Fe K a 6.7 keV emission-line-weighted abundances 
(middle column), and the mass- weighted temperature of the 
gas (right column), for five different redshifts from z ~ 1 to 
z 0.1. At each redshift, the cluster centre is defined as 
the centre of the most massive progenitor of the cluster at 
z = 0, and the distances from it are estimated using comov- 
ing coordinates. We can see that the ICM at z = is formed 
by the accretion and merger of small substructures that at 
early times (z ~ 1) are spread within an extended region 
of radius ~ 9/i _1 Mpc. These gaseous building blocks are 
already considerably contaminated by that time, with iron 
abundances that reach values of Fe/H ~ 0.3. 

From these figures, we can also infer interesting aspects 
about the dynamics of the contaminated gaseous clumps. 
Those that are at ~ 6/i -1 Mpc from the cluster centre at 
z ~ 0.5 needs to move at more than ~ lOOOkms -1 to get 
inside 0.5 h~ x Mpc by z — 0. A quantitative analysis of the 
dynamics of these gaseous clumps is presented in Figure 
which shows the dependence of the velocity of gas particles 
with their clustercentric distances and its evolution with red- 
shift. We identify the gas particles contained within spherical 
shells of 0.5 h~ x Mpc of thickness, centred at the most mas- 
sive progenitor of the cluster at z ~ 0.5, and estimate their 
mean velocities and clustercentric distances at lower red- 
shifts. In this way, we follow the evolution of these proper- 
ties as the cluster is being assembled. The different symbols 
denote the redshifts considered, and the lines connect the 
sets of gas particles defined according to their clustercentric 
distance at z ~ 0.5. There are many aspects to emphasize 
from this plot. At z ~ 0.5, gas particles have very high ve- 
locities, which range from - 1300 to - 2500km s" 1 . The 
highest velocities correspond to gas particles lying between 
~ 1 and ^ 2 Mpc. This is the same for redshifts z ~ 0.3 
and ~ 0.2, although the velocities acquired are smaller. If 
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Figure 5. Evolution of iron abundance and temperature of gas particles that lie within 500 kpc from the cluster centre at z = 0. 
The plots show the projection of mass- weighted iron abundance by number relative to hydrogen, Fe/H, respect to the solar value (left 
column), the corresponding Fe K a 6.7 keV emission-line- weighted abundance (middle column), and the mass- weighted temperature of 
the gas (right column). We give results for three redshifts: z ~ 1 (lookback time of 7.8 Gyr; first row), z ~ 0.5 (lookback time of 4.9 Gyr; 
second row), and z ~ 0.3 (lookback time of 3.6 Gyr; third row). At each refshift, the spatial coordinates are centred at the most massive 
progenitor of the cluster at z = and are expressed in comoving scales. 



we follow the evolution with redshift of gas particles lying 
more than 4/i _1 Mpc at z ~ 0.5, we see that they are accel- 
erated till they achieve the highest velocities at z ~ 0.3 and 
~ 0.2. This occurs at a distance corresponding to the virial 
radius of the main progenitor, which is of the order 1 h~ x 
Mpc at z rsj 0.5 and increases to ~ 1.5 h~ x Mpc at z = 0. 



Thus, despite shell crossing, it is evident that the dynamics 
of gas particles is characterised by a general behaviour in 
which they are accelerated till they reach the main cluster 
shock radius, at a distance of the order of its virial radius, 
and once they have been incorporated to it they converge 
to the centre with velocities of ~ 300kms _1 . This turbulent 
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Figure 6. Same as Fig. [5] but for redshifts: z ~ 0.2 (lookback time of 2.4 Gyr; first row), and z ~ 0.1 (lookback time of 1.2 Gyr; second 
row). 



velocity reached inside the virial radius and the infall that 
dominates t he gas outside i t are c onsistent with the results 
obtained bv ISunvaev et al.l (120031) based on the analy sis of 
cluster simulations at z — (|Norman Sz Brvanlll999r) . The 
complex dynamics of gas particles acquired during the hier- 
archical formation of the structure helps to explain the way 
in which the metallicity profiles develop, as we discuss later 
in this section. 

Figure |H1 shows the evolution of iron and oxygen abun- 
dance profiles, which are determined for gas particles lying 
within the innermost comoving 1 /i -1 Mpc at z — and at 
each of the redshifts considered in Figures |5] and |S] The first 
three rows show the shape of the profiles for the combined 
effect of both types of SNe, and for the separate contribu- 
tions of SNe la and CC, respectively. The plot in the last row 
gives the evolution of the O/Fe ratio profile resulting from 
the influence of the combined effect of both types of SNe. 
The way in which the iron and oxygen abundance profiles 
develop are quite similar, even though the former is mainly 
determined by the large amount of iron ejected by SNe la, 
while the latter is only determined by the SNe CC contri- 
bution. This fact leads to an almost flat O/Fe ratio profile 
at z — 0, with a slight increasing trend at large radii. The 
central abundances of both Fe and O, and the almost flat 



O/Fe ratio show that both types of SNe contribute in a sim- 
ilar fashion to every part of the cluster, that is, without a 
preference for SNe la to chemically enrich the inner regions 
of the cluster, and for SNe CC t o pollute its outskirts, a s 
inferred from some observations (jFinoguenov et al.1 2001). 
The main change in the profiles occurs in the central region 
(< 100 h~ x kpc) at late times by a progressive increase of 
the chemical abundances. 

The evolution of the radial abundance profiles of the 
main SNe la and CC products (iron and oxygen, respec- 
tively) from z <^> 1 to z — 0, can be explained by taking into 
account the history of metal ejection of the cluster galaxies 
as well as the dynamics of the gas associated to them while 
galaxies are being accreted onto the cluster. 

Figure shows the iron mass ejection rate as a function 
of redshift produced by SNe la (top panel) and SNe CC 
(bottom panel) that are contained in galaxies that lie within 
shells of different distances from the cluster centre at z — 0. 
Since we are focusing on a set of galaxies in the z — cluster, 
the contribution of SNe CC to the ICM chemical enrichment 
peaks at higher redshifts (z ~ 5 — 7) than the st ar formation 
rate o f the whole simulation (see figure 1 of ICiardi et alJ 
(2003) for a comparison of the redshift evolution of the star 
formation rate for field region simulations and for cluster 
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Figure 7. Mean velocities of gas particles contained within 
spherical shells of 0.5 h~ x Mpc of thickness, centred at the most 
massive progenitor of the cluster at z ~ 0.5, as a function of their 
mean clustercentric distances. This dependence is estimated for 
different redshifts as denoted by the symbols. The lines connect 
the sets of gas particles denned according to their clustercentric 
distance at z ~ 0.5. 



simulations). As we have already mentioned, the SNe CC 
rate closely follows the star formation rate since the lifetimes 
of the stars involved are quite short. Instead, the ejection 
rate of elements produced by SNe la reaches a maximum 
at z 3 — 5, due to the return time distribution resulting 
from the adopted model for SNe la explosions. The peaks in 
the mass ejection rates are followed by a strong decline at 
lower redshifts for both types of SNe, such that the ongoing 
chemical contamination is quite low at z = 0. This behaviour 
is more pronounced for those galaxies that are nearer to the 
cluster centre at the present time. 

Considering the accumulated mass of iron in the ICM at 
z = 0, we find that SNe la are responsible for ~ 75 per cent 
of the iron content of the ICM. Interpretation of observed 
profiles of abu ndance ratios extends this p ercentage up to 
~ 80 per cent JGastaldello fc Mo lendi 200%. However, both 
results are quite dependent of the assumed SNe la model. 
The point that we have to emphasize is that the relative 
contribution of SNe CC and SNe la to the iron content of 
the cluster is already one-third at z ~ 2, when the accu- 
mulated iron masses provided by both sources peack. The 
lack of metal contribution by SNe explosions at late times 
from galaxies that are near the cluster centre cannot explain 
the increase of abundances in the inner 100 h' 1 kpc since 
z ~ 0.1, suggesting that dynamical processes are playing an 
important role in the developement of abundance profiles. 

Figure 1101 shows the evolution of mass- weighted iron 
abundances of gas particles contained within l/i -1 Mpc 
from the cluster centre at redshifts z ~ 1, ~ 0.5, ~ 0.3, 
~ 0.2, ~ 0.1 and z — 0. We clearly see from these plots that 
highly enriched gas clumps are already present in the out- 




1.30L = 

0.1 1.0 
R [ h" 1 Mpc ] 

Figure 8. Evolution of the radial profiles of iron and oxygen 
abundances by number relative to hydrogen and of the profile 
of the O/Fe ratio, referred to the solar value. The profiles are 
determined for the gas contained within the innermost comov- 
ing l/i -1 Mpc at z ~ 1 (dash dot dot line), z ~ 0.5 (dash dot 
line), z ~ 0.3 (dashed line), z ~ 0.2 (dotted line), z ~ 0.1 (thin 
solid line) and z = (thick solid line). The first row shows the 
contributions of both types of SNe (CC and la) to the Fe and 
O abundances, while the second and third rows show the sepa- 
rate contributions of SNe la and CC, respectively. The last row 
presents the O/Fe abundance ratios when both types of SNe are 
considered. 

skirts of the cluster at z ~ 0.3 and that the iron abundance 
profile starts to settle down at z ~ 0.2. The contaminated 
gas clumps located near y ~ 350 h~ x kpc at ~ 0.2 and close 
to y ~ 200 /i -1 kpc at ~ 0.1 seem to converge to the cluster 
centre at z — considerably increasing the iron abundance 
in the inner 100 kpc. Thus, gas particles bound to the 
dark matter substructures of 'halo galaxies' (see definition in 
sect ion IZ2)> have been mainly enriched at early epochs, when 
the metal ejection rates from these galaxies were higher, and 
have then fallen to the inner regions of the cluster driven by 
the haloes of infalling galaxies. This dynamical effect is quite 
likely taking into account the analysis of gas dynamics based 
on figures El and A simple test of finding the position 
of the progenitor of the central galaxy at different redshifts 
reveals that these galaxies are at ~ 100 kpc at z ~ 0.5, 
reinforcing our conclusion based on gas dynamics. 

The scenario of chemical enrichment considered in 
our model indicates that gas dynamical effects play an 
important role in the developement of abundance pat- 
terns. However, there are other sources of chemical ele- 
ments not included in our model that may help to in- 
terpret e observational results, such as intracluster stel- 
lar populations and ram-pressure stripping. There is in- 
creasing evidence of the presence of i ntracluster stars 
based on observat ions of diffuse light (Feldmeier 2004; 
IZibetti e t_ajjj20 05l) and individual stars between cluster 
galaxies (lFeldmeierl20oilDurrell et all2002l : lGal-Yam et all 
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Figure 9. Iron mass ejection rate as a function of redshift pro- 
duced by SNe la (top panel) and SNe CC (bottom panel), con- 
tained in galaxies that lie within different distances from the clus- 
ter centre at z = 0: the dominant cluster galaxy (at the clus- 
ter centre, thick solid line), galaxies within the spherical region 
with radius 0.5 h~ x Mpc without the central galaxy (dashed line), 
galaxies within the spherical shells between 0.5 and l./i -1 Mpc 
(dash dot line), and between 1.0 and 1.5 h~ x Mpc (dash dot dot 
dot line). 

bOOflft . High resolution TV-Body /SPH simulations are start- 
ing to yield predictions of the spatial and k i nematic distri- 
butions for these stars { ffijlLmjm^ et alJ 
2004: ISommer-Larsen. Romeo PortinarilfeOO^T ) which hold 
important information about the assembly history of a 
galaxy cluster. These numerical works show that a fraction 
of at least ~ 10 per cent of unbound stars accumulate within 
the cluster as a result of stripping events and infall of galaxy 
groups that already contain unbound stars. The impact of 
intracluster stars in the chemical enrichment of the ICM has 
been explored by IZaritskv. Gonzalez &; Zabludoffl (|2004T) by 
using a simple model. They demonstrate that this stellar 
component makes a significant contribution to the iron con- 
tent of the ICM. The metal pollution due to ram-pressure 
stripp ing h as been recently i n vestig ated bv ISchindler et all 
( 2005|) and iDomainko et alJ (|2005T ) . These authors argue 
that the efficiency of this enrichment mechanism is higher 
than that of supernovae driven galactic winds, producing 
a centrally concentrated metal distribution in massive clus- 
ters. 



8 PROBING THE ICM DYNAMICS WITH 
METALS 

In the previous section we have showed how the complex 
dynamical evolution of intracluster gas influences the devel- 
opement of the abundance patterns of the ICM. We now 
discuss the potentiality of spectroscopic observations for de- 



tecting such strong gas bulk motions, sheding light on our 
understanding of gas dynamics during the cluster formation. 

The spatial and spectral resolution of X-ray telescope 
detectors onboard of present satellite missions allow the 
measurement of metallicity maps and of high resolution 
emission line spectra. The combination of these data con- 
tributes enormously to improving our knowledge about the 
ICM properties. For example, Fe K a 6.7 keV emission-line- 
weighted maps contain valuable information about thermo- 
dynamical and chemical properties of the intracluster gas, 
as discussed in Sectional We can also learn about the kine- 
matics of the intracluster gas and its spatial correlation with 
other properties from maps of emission-line-weighted line- 
of-sight (LOS) velocities and from Fe K a 6.7 keV emission 
line spectra generated by the gas along LOSs through the 
cluster. 

Figure ITT1 shows projections of Fe K a 6.7 keV emission- 
line- weighted LOS velocities of the gas particles onto the 
three orthogonal planes x-y, y-z and z-x. Velocities are rel- 
ative to the cluster centre, being colour-coded such that red 
and yellow areas correspond to gas moving towards the ob- 
server, and green and blue represent gas receding from it. 
Only material contained within a sphere of radius 1 h~ x Mpc 
centred on the dominant galaxy of the cluster has been in- 
cluded. These maps exhibit steep azimuthal and radial gra- 
dients, revealing a much more complex structure than those 
depicted by the density, temperature and iron abundance 
distributions, indicating the presence of large-scale turbulent 
motions generated by inhomogeneous infall. This is consis- 
tent with our results on the evolution of the dynamics of gas 
particles based on the analysis of Figure |7| The sharp edges 
denned by red and green colours near the cluster centre are 
contact discontinuities between gas that originated in differ- 
ent clumps. Note that this is not a projection effect of gas 
particles lying far from the cluster centre. These contact dis- 
continuities are also seen in the projection of particles well 
inside the cluster (< 50 kpc). 

The gradients visible in the velocity maps of Figure [H] 
can be as large as ~ 1000 kms -1 over distances of a few hun- 
dred kpc. They are of the same order of magnitude as those 
obtained from acc urate observations of the Centaurus clus- 
ter (Abell 3526) fchroke &: Bregmanl EoOlh . We note that 
this cluster is a very good candidate for this kind of stud- 
ies, being one of the closest X-ray bright clusters of galaxies 
with an optical redshift of 0.0104. Its core has been spatially 
resolve d by the Ch andra Observatory, using the ACIS-S de- 
tector (San ders fc Fabianl 2002V leading to temperature and 
abundance maps. However, there is not enough information 
in the spectrum to isolate spectral lines and perform accu- 
rate gas velocity measurements. This task is possible with 
the detector GIS on board of ASCA. Th e gas velocity distri- 
butio n determined with this instrument (|Duoke &; Bregmanl 
2001) shows a region associated with the subcluster Cen 45, 
at rsj 130 kpc from the main group Cen 30 (centred on 
the cD galaxy) , with a radial velocity higher than the rest of 
the cluster by ~ 1700 kms -1 . These radial velocity measure- 
ments obtained from X-ray observations support the opti- 
cally determined veloci ty difference bet ween Cen 45 and Cen 
30 of - 1300 kms -1 (IStein et al.lll997Tl . 

These gradients in radial velocities are observationally 
detected by the Doppler shifting of the metal lines used to 
perform the spectral fits, which are mainly driven by the Fe 
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Figure 10. Evolution of iron abundance of gas particles that lie within l/i -1 Mpc from the cluster centre at redshifts: z ~ 1 (upper 
left panel), z ~ 0.5 (upper middle panel), z ~ 0.3 (upper right panel), z ~ 0.2 (lower left panel), z ~ 0.1 (lower middle panel) and z = 
(lower right panel). The plots show the projection of mass- weighted iron abundance by number relative to hydrogen, Fe/H, respect to the 
solar value. At each refshift, the spatial coordinates are centred at the most massive progenitor of the cluster at z = and are expressed 
in comoving scales. 
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Figure 11. Projection of Fe K a 6.7 keV emission-line- weighted line-of-sight (LOS) velocities of gas particles contained within a sphere 
of radius 1 h~ 1 Mpc centred on the dominant cluster galaxy: x-y projection (left panel), y-z projection (middle panel), and z-x projection 
(right panel). 
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K a complex. However, a higher spectral resolution would al- 
low us to detect the different components in which a certain 
line is split due to the bulk motions of different magnitude 
along a line-of-sight through the cluster. Future X-ray mis- 
sions promise a substantial increase in energy resolution, 
paving the way for the construction of a full 3-D picture of 
the ICM dynamics. ISunvaev et alJ £2003) analyse this pos- 
sibility by calculating the synthetic spectra of the Fe K a 6.7 
keV emission line along sigh t lines through a simulated clus- 
ter (|Norman Brvanlll99d ). They show that this emission 
line is split into multiple components over a range of energy 
AE ~ ±15 eV, as a result of turbulent and bulk motions 
of the gas. They estimate that a spectral resolution of 4 
eV, similar to the one planned for future missions, would be 
enough for the detection of dynamical features in the ICM 
based on emission line spectra. 

Following the study made bv lSunvaev et all {2003), we 
analyse the way in which the special features of the gas mo- 
tion that become evident in the velocity maps of Figure 1111 
are reflected in the Fe K a 6.7 keV emission line spectra gen- 
erated by the gas along lines-of-sight through the cluster. 
The spectra are computed by dividing the sightlines in bins 
of width 20/i _1 kpc along directions perpendicular to the 
projection planes. The thermodynamical, kinematical and 
emission properties assigned to each bin are obtained from 
the 64 closest gas particles to the bin centres using the SPH 
smoothing technique. The energy of the iron line emitted 
from each bin along the sighline is Doppler shifted propor- 
tionally to the smoothed LOS velocity inferred from the gas 
around the bin. Combining these shifted energy values with 
the corresponding intensities of the line emission we obtain 
the iron emission line spectra along the LOS. On top of 
this, we consider the broadening suffered by each line due 
to thermal motions of the gas. 

Figures 1121 and 1131 show the resulting synthetic spec- 
tra for six LOSs along the z-axis through the x-y projection 
of the simulated cluster. Like in the construction of maps, 
the computation of the spectra only comprises gas particles 
within the virial radius of the cluster. Solid thick lines corre- 
spond to spectra that are affected by both Doppler shifting 
and thermal broadening. For comparison, we also include 
in each panel the iron emission line affected only by ther- 
mal broadening (dash-dotted lines). The maximum values 
of these last ones are used to normalise both kinds of iron 
emission lines, with and without Doppler shifting. The ra- 
tios of these maxima with respect to the one corresponding 
to the cluster centre, referred to as / max //™ ax 5 are a i so given 
in each of the panels that show the spectra, together with 
the location where the LOSs intercept the ICM. Below these 
panels, both FiguresH21andlT31displav the distribution of gas 
properties along the corresponding sightlines: emission-line 
weighted LOS velocities used to generate the spectra, and 
mass- weighted hydrogen density, temperature, iron abun- 
dance and normalised emissivity of the K Q 6.7 keV emission 
line. These distributions allow us to see how gas compo- 
nents characterised by different thermodynamical and chem- 
ical properties move at different velocities, as well as which 
parts of the line of sight contribute most to the spectral 
lines. 

An inspection of these distributions indicates that the 
gas reaches high iron abundances at high density and low 
temperature regions, as it is also shown by the maps of 



Figure El This situation is especially clear in the cluster 
core and in a subclump moving at high speed located at 
~ 500 /i -1 kpc from the cluster centre, as it is evident from 
the first and second columns of Figure The Doppler 
shifting and split of the emission line at the cluster core 
are produced by the large-scale turbulent motions of the 
gas, which produces a LOS velocity gradient from ~ 200 to 
—400 km s _1 around the middle of the sightline. This motion 
is reflected in the spectra because the intensity of the iron 
emission line is larger than 10 per cent the central value 
within a radius of ~ 100 /i -1 kpc from the cluster centre. 
However, this intensity decays quite abruptly at larger radii 
and becomes undetectable (see also the mass- weighted emis- 
sion map in Figure 0. For this reason, the subclump mov- 
ing at — 800kms _1 from the cluster centre does not pro- 
duce any signature in the first two Doppler shifted spectra 
shown. If the emissivity were higher at the location of this 
substructure, a peak would appear at ~ 6, 685eV. The sec- 
ond and third spectra lines in figure are produced by 
sightlines that intercept the cluster at a distance of 100 and 
360 /i -1 kpc from its centre, respectively. Note how the ratio 
jmax^jmax decreases as we get further from the cluster core. 
The second spectrum shows a clear Doppler shifting and 
broadening of the line, while the third one is only slightly al- 
tered with respect to the reference broadened line, in agree- 
ment with the LOS velocity distributions within the inner 
- 100 h' 1 kpc. 

Figure H31 presents three more synthetic spectra gener- 
ated from lines of sight that intercept the cluster through 
areas with large offset velocities near its centre (see fig- 
urelTTl. The first one shows a split of the emission line in two 
components arising from the central velocity gradient from 
~ —200 to 100 km s -1 . The line in the second spectrum is 
widely Doppler shifted as a result of high bulk motions of 
~ 500kms _1 in the central parts of the sightline. In the 
third emission line spectrum, a second component becomes 
apparent as a result of gas moving at ~ — 500 km s -1 . As we 
noted in the previous set of spectra in Figure the strong 
gas motions of ~ — 500kms _1 and ~ — 700kms _1 occurring 
at rsj 500 h~ x Kpc from the centre of the lines of sight con- 
sidered in the second and third columns are not imprinted in 
the corresponding spectra because of the low emissivity that 
characterises this region. The shifts and splits observed in 
all these spectra are of the order of ~ 10 — 15 eV. This effect 
is too small to be resolved by present instruments, since, for 
example, the spectral resolution of the EPIC pn detector for 
the Fe K a 6.7 keV emission line is 150 eV. 

From this analysis we can conclude that only gas mo- 
tions that produce absolute values of LOS velocities higher 
than 400kms _1 around the cluster core, where the iron 
line emissivity is larger than 10 per cent the central value, 
could become observable by future X-ray missions through 
the shifting and split of the line in two components due 
to Doppler effects. Lower energy lines like the Fe L-blend 
centred around 1 keV could be also used for the present 
purpose. In principle, they should show much more struc- 
ture since they are strongly emitted in lower temperature 
gas. However, gas with this characteristic lies far from the 
centre of our hot cluster and the low line emissivity in these 
areas prevent us from getting more information from these 
lines than that provided by the Fe K Q line at 6.7 keV. These 



Metal Enrichment of the I CM 19 



I. [km s" 1 ] LOS vel. [km s" 1 ] LOS vel. [km s" 1 ] 

C) 500 -500 u 500 -500 u 




-1000 -500 500 1000 -1000 -500 500 1000 -1000 -500 500 1000 

Position along LOS [h 1 kpc ] Position along LOS [h 1 kpc ] Position along LOS [h 1 kpc ] 



Figure 12. Fe K a 6.7 keV emission line spectra along three lines of sight (LOSs) on the z-axis through the simulated cluster, which is 
projected onto the x-y plane, and the distribution along the corresponding LOS of kinematical, thermodynamical, chemical and emissivity 
properties of the gas that contributes to the spectra. The upper panels show the effect of thermal broadening on the Fe emission line 
(thin dash-dotted curve) and the composition of thermal broadening and Doppler shift due to gas bulk motion with respect to the cluster 
centre (thick solid line). The maximum values of the broadened lines are used to normalise both kinds of iron emission lines, with and 
without Doppler shift. These panels indicate the ratio of this maximum with respect to the one corresponding to the cluster centre, 
7 max //™ ax , and the x-y position where the LOS intercepts the ICM. The five plots below each spectra show the distributions of gas 
properties along the corresponding sightlines: emission-line weighted LOS velocities used to generate the spectra, and mass-weighted 
hydrogen density tih, temperature T, iron abundance Fe/H, and normalised emissivity of the K a 6.7 keV emission line 7// max . Velocities 
are expressed in kms -1 , wh in cm -3 , T in K, and iron abundances are referred to the solar value. The positions along the LOSs are 
given with respect to the cluster centre. 
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Figure 13. Same as figure IT2l for other set of LOSs that intercept the ICM through areas with large offset velocities present near the 
cluster centre as depicted by Figure HT1 



lines are probably much more usefull in a lower temperature 
cluster. 

We would like to note that the hydro- dynamical sim- 
ulation considered by ISunvaev et alJ ()2003h ignores radia- 
tive cooling, which affects the thermal structure in cluster 
cores, and they assume constant iron abundances to calcu- 
late the iron emission line spectra. Our simulation suffers 
from the first drawback as well; however, the critical advan- 
tage of our hybrid model over approaches like the one used 



by ISunvaev et al.l (|200rf ) is that the chemical abundances 
of gas particles are dynamically and consistently generated 
from the stars in the galaxies. 

Taking into account the information supplied by LOS 
velocity maps and by spectra along selected sight lines, we 
can connect the global features of bulk motions in the gas 
with the multiple components in the spectra that originate 
as a result of these motions. Note that the latter feature 
can be more easily observationally detected and quantified. 
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Thus, our results can constitute valuable help for the inter- 
pretation of observational data. 



9 SUMMARY AND CONCLUSIONS 

We have presented an hybrid model for the chemical enrich- 
ment of the intracluster gas that allows a detailed analy- 
sis of the spatial distribution of chemical, thermodynamical 
and kinematical properties of the ICM. Our model combines 
a cosmological non-radiative hydrodynamical iV-Body/SPH 
simulation of a cluster of galaxies, and a semi-analytic model 
of galaxy formation. The spatial distribution of metals in 
the ICM can be accounted for by the chemical enrichment 
of the gas particles in the underlying AT-Body/SPH simu- 
lation. This enrichment results from the metals ejected by 
the galaxy population, which are generated by applying the 
semi-analytic model to the the dark matter haloes detected 
in the outputs of the simulation. This link between semi- 
analytic model results and the chemical enrichment of gas 
particles leads to a consistent chemical enrichment scheme 
for the intracluster gas, being the key ingredient to carry 
out the present project. 

The parameters that characterise our model are chosen 
according to several observational constraints. The proper 
circulation of metals among the different baryonic compo- 
nents leads to a spatial distribution of chemical elements 
in the intracluster gas that closely resemble observed abun- 
dance patterns, as can be seen from the comparison of the ra- 
dial abundance profiles with data from lTamura et ail (|2004) 
(Figure 0. The corresponding mean iron abundance of the 
intracluster medium at z = is [Fe/H] ~ 0.28. It has been 
reached with the contribution of SNe la that provide ~ 75 
per cent of the total ICM iron content at z = 0. While this 
percentage is quite dependent on the parameters that char- 
acterise the chemical model implemented, what reflects the 
circulation of metals is the fraction of elements from differ- 
ent sources that end up in the ICM as opposed to in stars. 
The fraction of iron originated in SNe la that is contained 
in the ICM is ~ 0.78 while that locked in the stars of cluster 
galaxies is ~ 0.21. These fractions are slightly different for 
the iron ejected by SNe II, being ~ 0.72 and ~ 0.28 for hot 
gas and stars, respectively. The increment of the proportion 
of SNe II products in stars reflects the relative time-scales 
for star formation and metal pollution by different types of 
super novae. 

After showing that our model reproduces reasonably 
well the distribution of metals in the ICM at z = 0, we anal- 
yse the evolution of chemical and dynamical properties of 
the intracluster gas in order to understand the physical pro- 
cesses involved in the determination of its final abundance 
patterns. The aim of this study is also to provide hints that 
helps in gaining information about the characteristics of gas 
bulk motions from spectroscopic observations. We analyse 
the rich information provided by our hybrid model in three 
ways, 

(i) by constructing radial abundance profiles of the ICM, 

(ii) by building projected mass- weighted and Fe K a 6.7 
keV emission-line-weighted maps of chemical, thermody- 
namical and kinematical properties of the diffuse gas, 

(iii) and by calculating synthetic spectra of the Fe K a 



6.7 keV emission line along sight lines through the simulated 
cluster. 

The spatial distribution of chemical elements leads to 
radial abundance profiles with negative gradients, irrespec- 
tive of the type of source from where they originated. The 
central enhancement of the Fe and O abundance profiles and 
the flat profile of the O/Fe ratio support the notion that 
the intracluster gas is polluted in the same way by both 
types of SNe. From the analysis of the evolution of these 
radial abundance profiles and the projected mass- weighted 
iron abundance of gas particles at different redshifts we find 
that a high level of enrichment of the ICM is already reached 
at z rsj 1. Despite the delay time for the ejection of metals 
inherent to SNe la, the contribution of both types of SNe to 
the ICM chemical enrichment peaks at z ~ 2. 

The lack of evolution of mean iron abundances up to 
z ~ 1 is in agreem ent with observations of distant clusters 
dTozzi et alJ l2003V However, the final abundance profile of 
the cluster gets steeper within the inner ~ lOO/i" 1 Kpc dur- 
ing the last Gyr (see the increments between z ~ 0.1 and 
z = in Figure This cannot be accounted for the metal 
production of the cluster galaxies. The peaks in the rates 
of the ejected iron mass from these galaxies are followed 
by a strong decline towards lower redshifts for both types 
of SNe, such that the ongoing chemical contamination at 
z = is quite low; this behaviour is more pronounced for 
galaxies that are nearer to the cluster centre at the present 
time. The combination of all these results indicates that the 
enhancement of central abundances produced at low red- 
shift cannot be explained by metals ejected during this time 
interval from galaxies that lie close to the cluster centre. 
Indeed, galaxies within a sphere of radius 500 kpc cen- 
tred on the dominant cluster galaxy do not contribute to the 
metal enrichment at recent epochs. Instead, a scenario ap- 
pears favoured in which gas particles are primarily enriched 
at high redshifts, and are subsequently driven to the cluster 
centre by bulk motions in the intracluster gas. Hence, the 
main cause that contribute to this central abundances en- 
hancement is the dynamical evolution of gas particles during 
the cluster formation. 

Gas particles develop very high velocities ranging from 
~ 1300 to - 2500 krns" 1 when they are farther than 
~ 1 h~ x Mpc from the cluster centre at z > 0.2. Once 
these high velocity gas clumps cross this typical cluster 
shock radius, which correspond approximately to the virial 
radius of the most massive progenitor of the cluster at 
these redshifts, they slow down ending up within the inner 
~ 500/i _1 Kpc of the cluster at z = 0. The inhomogeneous 
gas infall produces contact discontinuities that are mani- 
fested as sharp edges in the emission-line-weighted line-of- 
sight velocity maps. These maps provide useful information 
about the gas dynamics in the ICM and its dependence on 
spatial position within the cluster. They reveal a quite com- 
plex structure of the gas dynamics, including the presence 
of strong radial and azimuthal gradients with values as large 
as rsj lOOOkms -1 . Gradients of this magnitude are also seen 
in velocity distributions along sightlines through the clus- 
ter centre. Doppler shifting and broadening suffered by the 
Fe K a 6.7 keV emission line along sightlines could be used 
to probe these gas bulk motions when they are produced 
within an area characterised by high iron line emissivity. 
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Velocities of ~ 400 km s -1 , typically found within a radius 
of 100 h~ x kpc from the cluster core, produce Doppler shifts 
of the line of the order of ~ 10—15 eV. The line can be split 
in multiple components depending on the structure of the 
velocity field. We find that the higher velocities reached by 
the intracluster gas beyond 500 hT 1 kpc, or those aquired by 
infalling subclumps located at such large distances from the 
cluster centre do not produce a detectable signature in the 
spectral lines because of the rather low iron line emissivity 
in the outskirts of the cluster. 

The combination of all these results strongly supports 
the fact that abundance patterns of the ICM are the result 
of gas mixing because of hydrodynamical processes during 
cluster formation, being a considerable amount of the gas 
mainly enriched at high redshifts, before the cluster have 
time to virialise. The central enhancement of Fe and O abun- 
dance profiles might be considered as produced by a partic- 
ular situation taking place in our simulation, in which two 
big highly enriched gas subclumps converge at the cluster 
centre at the present epoch. However, this kind of behavi- 
uor is consistent with the hierarchical formation scenario, in 
which accretion and merging of substructures constitute the 
main processes driving the formation and evolution of clus- 
ters of galaxies. Thus, the specific shape of metallicity pro- 
files might well be a consequence of the combination of two 
time-scales, one related to the metal production of galaxy 
clusters and the other to the metal mixing. In this last case, 
both the effect of gas dynamics, galaxy mergers and the effi- 
ciency in the diffusion of metals due to energetic and chemi- 
cal feedback associated with the star formation process play 
an important role. Other enrichment processes not taken 
into account in our model, such as metal pollution by in- 
tracluster stars and ram-pressure stripping, may contribute 
significantly to recover the observed abundance patterns re- 
inforcing the effect of gas dynamics discussed here. 

Our analysis techniques in this study yield a 3-D pic- 
ture of the chemical and dynamical properties of the ICM. 
The information about the latter can be extracted from the 
metal content of the intracluster gas. Our results show that 
gas bulk motions are imprinted in the shape of X-ray lines of 
heavy ions. Since thermal broadening is small, these features 
could be detected with the high-resolution spectrographs of 
future X-ray missions (CONSTELLATION-X and XEUS). 
Such observations should be very important in determining 
the physical processes involved in the formation and evolu- 
tion of galaxy clusters, contributing with clues that allow us 
to verify the enrichment scenario inferred from our model in 
which gas dynamics has a crucial effect. 

Finally, even though all the processes considered in our 
hybrid model are tuned to satisfy numerous observational 
constraints, the implied global level of iron enrichment we 
find is on the low side of the generally accepted observed 
range, which is around one-third the solar value. This may 
in part be due to the uncertainties in the stellar yields and to 
our poor knowledge of the strength and nature of the feed- 
back processes involved, but could also point to a systematic 
problem with the observational estimates. Our model so far 
only accounts for the effect of SNe driven outflows, but it 
should be possible to include other sources of feedback, like 
those from AGN, as well. This is left for future work, and 
may perhaps provide a solution to some of the puzzles still 
present in understanding the observed cluster metallicities. 
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